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Abstract. The feasibility of a measurement of the electric and magnetic nucleon form factors at B-meson 
factories through the radiative return is studied. Angular distributions allow a separation of the contri- 
butions from the two form factors. The distributions are presented for the laboratory and the hadronic 
rest frame, and the advantages of different coordinate systems are investigated. It is demonstrated that 
values up to 8 or even 9 GeV^ are within reach. The Monte Carlo event generator PHOKHARA is extended 
to nucleon final states, and results are presented which include Next-to-Leading Order radiative corrections 
from initial-state radiation. The impact of angular cuts on rates and distributions is investigated and the 
relative importance of radiative corrections is analysed. 



1 Introduction 



The importance of measuring the nucleon form factor has 
been repeatedly emphasized in the literature (e.g. 
EI and refs therein). Recent experiments at Jefferson Lab 
have explored the ratio of electric and magnetic form fac- 
tors in the space- like region |3S1 up to 5.6 GeV^, using 
the recoil polarization method; their results show disagree- 
ment with data obtained with the Rosenbluth method 
(For an extensive review see [H). For recent suggestions 
of a possible explanation of this discrepancy through cor- 
rections from two-photon exchange see e.g. |ZIIHI-)- Data in 
the time-like region based on electron-positron annihila- 
tion into proton-antiproton (and neutron-antineutron) fi- 
nal states, and the inverse reactions, only extend to 6 GeV^ 
and 14 GeV^ respectively. These latter measurements ex- 
hibit fairly large errors and do not provide a separation 
of the contributions from the two form factors. To explore 



the validity of the different model predictions, improved 
measurements over a wide kinematic range are desirable. 

In the present paper the potential of the radiative re- 
turn at B-meson factories will be explored. As shown be- 
low, these measurements will cover the region from thresh- 
old up to 8 or even 9 GeV^ with sufficient counting rates. 
This region is also accessible at the Beijing storage ring 
and is complementary to the one at Cornell with energies 
above the J/4' resonance. 

Not surprisingly, many aspects of the radiative return 



are quite similar to those of the reaction e^e 



pp. In 
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particular it is again the square of the electric and mag- 
netic form factors, which can be determined separately 
through the analysis of angular distributions. 

Radiative corrections are indispensable for a precise 
interpretation of the data. Furthermore, given the limited 
acceptance and the asymmetric kinematic configuration at 
present i?-meson factories, a Monte Carlo event generator 
is required to demonstrate the feasibility of the measure- 
ment. The present analysis is based on an extension of 
the generator PHOKHARA 3.0 which was originally 
constructed to simulate two-pion and two-muon JO], and 
later four-pion ^11A2, . production through the radiative 
return |13lll4j and which includes next-to-leading order 
(NLO) radiative corrections |15II16| . The new version of 
PHOKHARA (PHOKHARA 4.0) will include, besides of 
the nucleon-antinucleon final states, final-state radiative 
corrections (FSR) to production at NLO |17j . 

For the present case we consider photon emission from 
the initial state only (ISR). As discussed below in more 
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detail, final-state radiation is completely negligible if the 
reaction is considered at a centre of mass (cms) collider 
energy of 10.5 GeV. 

The choice of form factors made by the program is di- 
rectly taken from [Tlll8| . It is, however, set up in a modular 
form such that the two form factors can be easily modified 
and replaced by the user, if required. 

The program also allows a more detailed study of an- 
gular distributions, including the effects of cuts on photon 
and proton acceptance. Indeed we will find that these dis- 
tributions are quite sensitive to the ratio of the form fac- 
tors, and the extraction of this ratio seems to be feasible 
over a fairly large kinematic range. 



2 Electric and magnetic form factors and the 
radiative return 

The matrix element of the electromagnetic current gov- 
erning nucleon-antinucleon production is conventionally 
written in terms of the Dirac (F^) and Pauli (F^) form 
factors 



Ja = -^e ■ u{q,) ( F^ {Q^)^, - [7,, 1 



where N stands for proton or neutron. The nucleon and 
antinucleon momenta are denoted by qi and q2 respecti- 
vely, and Q = qi + q2- 

Cross sections and distributions for e+e^ NN are 
conveniently expressed in terms of the electric and mag- 
netic Sachs form factors (191: 

G^, = F,^ + F,^, G%=F^ + rF^, (2) 

with T = Q^/Amf^, which will also lead to compact for- 
mulae for the radiative return. Both proton and neutron 
form factors can be decomposed into isoscalar and isovec- 
tor contributions (see e.g. (18;): 



(3) 



Let us start with a qualitative discussion of form factor 
measurements through the radiative return, based on the 
leading order process 



e+(pi) + e-(p2) ^ N{qi) + N{q2) + 7(fc) 



(4) 



From the simple analytical results it will be straightfor- 
ward to evaluate production rates, to understand the qual- 
itative aspects of angular distributions, and to develop 
strategies for the separation of the electric and the mag- 
netic form factor. 

The differential cross section for reaction Q (with ISR 
only) can be written as 



da 



1 



L^,^Hi'''d^2{pi +P2;Q, k)d^2{Q\ qi,q2 



dQ^ 
' 2tt 



(5) 

where L^j, and H^'^ are the leptonic and hadronic ten- 
sors respectively. For notation, definitions and an explicit 



form of the leptonic tensor, see for instance |12lll()| . The 
hadronic tensor is given by 



F[i_iu — 2|G]J^p((5^(5iy — g/iuQ 



-|G^n<ZM9-, (6) 



where q = (92 ^ From the explicit form of the 

hadronic tensor it becomes apparent that the relative phase 
between and GjJ^ cannot be measured in this experi- 
ment (i.e. as long as the spin of the nucleon remains un- 
measured). This is independent of the detailed form of 
the leptonic tensor. In particular phases from higher or- 
der ISR or (longitudinal or transversal) beam polarization 
affect only the leptonic tensor and thus do not alter this 
conclusion. 

Integrating over the whole range of nucleon angles and 
the photon azimuthal angle, and restricting the photon 
polar angle within ^ < 6i;°'" < 9^ < tt - 6^'"^, the dif- 
ferential cross section factorizes into the cross section for 
electron-positron annihilation into hadrons and a flux fac- 
tor that depends on s and only ;13) : 



da Aa'- 



(I) dg2 3sQ' 



1 ^ , 

1 c„ 



where c„i = cosf?™™, rrie denotes the electron mass, 



(7) 



i?(g' 



(7(e^ 



NN) 



f^point 



f3N_ 

2 



^IG^A' 



with cTpoint the lowest order muonic cross section and f3 



(8) 



N — 

1 - 4m%/Q\ 

After integration over the whole range of nucleon an- 
gles, the separation of electric and magnetic form factors 
is no longer feasible. However, the modulus of their ratio 
can be determined from the study of properly chosen an- 
gular distributions. The fully differential cross section is 
essentially given by 



32s 



G 



N |2 



N\2 



(9) 



)Vyi 



2 IG 



N |2 
Ml 



IG 



^12 



1 

Z/2 
1 1 

v\ y2 



ipi ■ qf + {P2 ■ qf 



s(s-g2) 



where yi,2 = ^-57- (1 T cosOj). 

The separation of the form factors will therefore have 
to rely on the peculiar dependence of the differential cross 
section on the nucleon and antinucleon momenta sepa- 
rately. 

The close connection between this analysis and the one 
based on the reaction e+e~ NN becomes even more ap- 
parent once the result is expressed in terms of the polar 
and azimuthal angles and tp, which characterize the nu- 
cleon direction in the rest frame of the hadronic system, 
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with the z-axis opposite to the direction of the real pho- 
ton momentum and the y-axis in the plane spanned by 
the beam and the photon direction (see Appendix A): 



(10) 



G 



N |2 
J\/l 



\G 



N\2 



4Q2 



1 

yi 



1 

y2 



Q2 \V r'^';(s-Q2 

X I {(3"f cos OY + ["i cos 0^ cos 9 — sin 9^ sin 9 sin (p)^ 



2 \G 



N |2 
Ml 



1 



1 

y2 



where 7 = (s + Q'^)/2y^ and ^ = (s - Q'')/{s + Q^) 
characterize the boost from the laboratory to the hadronic 
rest frame. In the limit <C s, this can be approximated 

by 



(47ra)3 (1 + cos2 6*^) 

Q2 (l_cos2 6'^) 



(11) 



x4(|Gf,p(l + 



cos 



An alternative choice of the coordinate system, which 
is obtained through an additional rotation around the x- 
axis in the hadronic rest frame, leads to a diagonal form of 
the leptonic tensor. In this frame the angular distribution 
of the baryons simplifies even further. This formulation is 
discussed in detail in the appendix. For small Q'/s the 
two frames nearly coincide. 

It is instructive to compare the angular distribution in 
Eq. 1)11(1 with the angular distribution from e+e~ NN: 



rate for baryon pair production from J/ip decays will be 
significantly enhanced, all considerations given above for 
continuum production do apply, and a precise measure- 
ment of the branching ratio of J/^ into baryon pairs and 
the relative strength of magnetic and electric coupling is 
within reach. 



3 Nucleon form factors and their 
implementation in PHOKHARA 

The parametrization of the form factors used in this paper 
and listed below in detail is adopted from Ref. with 
the analytical continuation to the time-like region taken 
from and The isoscalar and isovector components 
of Dirac and Pauli form factors are given by: 



— 



pv 



(0.120 

2 

g{Q'' 



(1 - -Pp-Tp 
' 3.706 • To 



with 



T = 



ml + (Q2 - ^ml)rpa{Q'^)lm^ 



(13) 

(14) 
(15) 
(16) 

(17) 



da 



N 



4Q2 



r<N \2 



(1 



-IG 



Af |2 



sm 



(12) 

The close relation between Eqs. (|11|) and (|12|l is clearly 
visible. 

As is already evident from the hadronic tensor in Eq. © , 
the phases of the form factors are, however, not accessible 
by this method. Their determination would require the 
measurement of the polarization of the nucleons in the 
final state PDI2T] . PHOKHARA can easily be extended 
to describe this additional degree of freedom. This option 
is of particular interest for AA production, with its self- 
analysing decay mode. 

As emphasized above, only ISR is included in the pre- 
sent version of the program. FSR with one photon only 
emitted from the hadronic system is described by the am- 
plitude for e^e~ ^ 7* ^ NNj, with an invariant mass of 
the virtual photon of 10.5 GeV. This option is, however, 
strongly suppressed by the proton form factor. In fact, 
the arguments are quite similar to those for pion pair pro- 
duction, where FSR in leading order and at 10.5 GeV is 
completely negligible. 

As will be discussed below, the event rate drops rapidly 
with increasing and only few events will be observed 
for invariant hadron masses around 3 GeV. However, the 



a(Q2 



and 



Tuj.fj) — 



1/2 



■log 



Q' 



, g{Q^) = 



l + \/]~a?2 



1 

(l-7e^»g2)2 



(18) 



2m^ 



The values of the parameters (dimensionful quantities 
in units of GeV) are [3p = 0.672, (5^ = 1.102, f3^ = 0.112, 
7710 = 1.019, rup = 0.765, m„ = 0.784, = -0.052, 
Fp = 0.112, 7 = 0.25. The angle 9 in Eq. lO is set to 
6' = 7r/4or6' = 197r/64 as recommended in J]. More recent 
fits to the data in the space-like region are listed in |22ll23j . 
However, in view of the large uncertainty of the data the 
discrimination between the different parametrizations is 
not yet possible. It should be emphasized that the original 
fit of ^18, agrees well with the ratio of the form factors 
measured at Jefferson Lab as shown in Fig.^ In the 
time-like region, the ratio of the form factors is predicted 
to increase with (Fig.|21), in contrast to the behaviour 
in the space-like region. 

Data giving information about the form factors in the 
time-like region originate from the reactions e^e~ — > pp^ 
pp —> e~^e~ and e~^e~ — > nn and are shown in Figs. OHHl 



4 



Henryk Czyz et alNucleon form factors, B-meson factories and the radiative return 




1 2 3 4 5 

-t (GeV^) 



Fig. 1. The ratio of proton form factors in the space-hke region 
as predicted by the model (sohd hne) and data from JLab 
|4I5| . fip is the proton magnetic moment. 




3.5 4 45 5 5.5 6 6.5 7 7.5 8 



s (GeV2) 

Fig. 2. Modulus of the ratio of proton form factors in the 
time- like region as predicted by the model |18|. 

Only measurements of the total cross section are available. 
Form factors, if available at all, are extracted under the 
assumption = /ZpG^, and no true separation of the 
electric and magnetic form factors has been performed, a 
consequence of the low event rates. Within the large exper- 
imental uncertainties the data are in reasonable agreement 
with the model predictions. The poor experimental knowl- 
edge can be substantially improved by using the radiative 
return |18II14| at fi-meson factories. Even the separation 
of electric and magnetic form factors may be envisaged in 
order to resolve the present discrepancies (H] in the mea- 
surement of the proton form factors. 

The form factors defined in Eqs. H13(l ~ (|16|l were im- 
plemented in the event generator PHOKHARA fTTllMTn] . 
providing a powerful tool for the experimental and theo- 
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Fig. 3. Comparison of the measured |Mli^Sli^6li^'7li^^lMl^ft| 

e^e~ —> pp cross section with the model from Ref. (181 . Pre- 
dictions are given for two different values (7r/4 - curve A and 
197r/64 - curve B) of the parameter 9 (see Eq. 1191 '). 
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Fig. 4. Comparison of the measured '31"32*33'34"35' pp 

e^e~ cross section with the model from Ref. 18 . Predictions 
are given for two different values (7r/4 - curve A and 197r/64 
- curve B) of the parameter 9 (see Eq. 11911 ). The curves were 
obtained without imposing any angular cuts. Theoretical re- 
sults with cuts appropriate for each experiment are shown as 
crosses. 



retical analysis of the processes e^e^ ~^ pp^ and nnj at 
NLO. Real and virtual FSR contributions are not (yet) 
implemented. They might, however, become important in 
conjunction with hard ISR (lowering the effective Q^) 
for close to the pp threshold, where the final-state 
Coulomb interaction is sizeable. 

The technical correctness of the implementation was 
checked first by confronting the distribution of events 
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Fig. 5. Comparison of the measured |24I26| e+e" nn cross for two different sets of cuts, 
section with the model from Ref. jl8j . Predictions are given for 
two different values (7r/4 - curve A and 197r/64 - curve B) of 
the parameter 9 (see Eq. 1191 '). 



generated by PHOKHARA 4.0 at LO with the correspond- 
ing analytical prediction in Eq. j?)). The test was success- 
fully performed with an agreement better than 0.1% for 
< 8 GeV^ (the region were the cross section is well ac- 
cessible to experiment). At NLO the standard tests for the 
independence of the cross section on the 'soft-hard' sep- 
aration parameter were performed again with a precision 
better than 0.1%. They guarantee that the implementa- 
tion of the form factors in the one-photon and two-photon 
parts of the generator is consistent and, together with the 
LO test, support the correctness of the whole implemen- 
tation. 



4 Form factor measurements at S-meson 
factories 
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Fig. 7. Differential cross section of the process e^e~ — > ^717(7) 
for two different sets of cuts. 



The differential cross section as a function of Q^, as ob- 
served through the radiative return, is shown in Figs. El 
and[3 The upper curve represents the cross section after 
integration over all proton angles and over photon angles 
down to 5° w.r.t. the beam direction. The lower curves 
are valid for cuts corresponding to the BaBar acceptance 
region transformed to the e+e" cms. The cross sections 
corresponding to the lower curves, integrated over < 
10 GeV^, amount to 59.3 fb for protons and to 125 fb 
for neutrons. With a luminosity of over 130 fb"-'^ accumu- 
lated at i?-factories already now, the event rate may be 
large enough to extract and separate \G^.j\ and \G^\ as 
proposed above. 

In Fig. |S1 we show the angular distributions in 0q (the 
polar angle of the vector q) for different ranges of Q^. 
To explore the sensitivity to the form factors individually 
we compare the differential cross sections obtained for the 
model described above with the distribution obtained un- 



der the assumption G^j ~ ^J-pG^, subject, however, to 
the constraint that a{e^e^ pp) (see Eq. remains 
unaffected. Without that constraint the predicted cross 
section would be up to ten times larger, in disagreement 
with existing data. The distributions are markedly differ- 
ent and the discrimination between different assumptions 
about the ratio G^z/G^ seems feasible. 

The angular distribution of the baryons, if defined in 
the laboratory or cms frame, is strongly affected by the 
boost. The difference between different choices for the 
form-factor ratios is expected to be more pronounced for 
the proton angular distribution in the hadronic rest frame, 
with the z-axis aligned with the direction of the photon 
and the y-axis in the plane spanned by the beam and the 
photon directions (see Eqs. (|10|1 and 1)11(1 above). In this 
case one may study alternatively the distributions with 
respect to 9 and ip, the azimuthal and polar angles of the 
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Fig. 9. Angular distributions in the polar angle of the pro- 
ton for three different ranges of Q^, with and without the con- 
straint ff^j — lJ.pG^ (see text for a more detailed explanation). 
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Fig. 11. Comparison between LO and NLO predictions for the 
differential cross section of the reaction e^e~ — > p|)7 for two 
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proton directions. For <C s, the case relevant to the 
present discussion, the ip dependence is unimportant, and 
the pronounced dependence of the cos distribution on 
the choice of the form factors is clearly visible. The differ- 
ence between the two model assumptions is significantly 
more pronounced, and the sensitivity to the form-factor 
ratio improves (Fig. ISJ . 

A quantitative estimate of the precision of such a cross 
section measurement for proton production can be ob- 
tained from Fig. 6. As is evident from Eq. (7), the relative 
error in the distribution da/dQ^ for a specific Q^-interval 
is a direct measure of the relative error in R{Q^), averaged 
over the same interval. Taking, for example, the angular 
cuts for photons and protons which roughly correspond 
to the detector acceptance, an integrated luminosity of 
100 fb~^ and a bin size of 0.5 GeV^, one expects about 
2500 / 600 / 250 events around 4/5/6 GeV^ and a 
resulting statistical precision of 2 / 4 / 6 %. For 1000 fb~^ 
even 8 GeV^ could be reached, again with roughly 6% sta- 
tistical precision. For neutrons the rates are even higher. 



The sensitivity to the ratio |Gm/(mG'£;)| can best be 
estimated from Fig. 9. Let us discuss the interval between 
4.5 to 5 GeV^. We expect in total around 1000 events, and 
the discrimination between the two options shown in the 
Figure should be straightforward. 

The importance of the radiative corrections can be de- 
duced from Figs. 10 and 11. Even if the integrated cross 
sections are similar in LO and NLO (a difference below 
0.5% is observed in both Figs. 10 and 11), the radiative 
corrections do lead to distortions of the distributions by 
1-3%; furthermore, they depend on the details of the cuts 
and the criteria for event selection. Hence, the use of the 
NLO generator is highly recommended. 

5 Summary 

The radiative return at B-meson factories is well suited 
for measurements of the nucleon form factors over a wide 
kinematic range. It should be emphasized that these mea- 
surements of the proton and neutron form factors can be 
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obtained from the data sample taken in standard runs 
and close to the T-resonance. Valid for the time-like re- 
gion, they would complement the precise results for the 
space-like region from JLab. Close to the threshold the 
statistical precision would be at the per cent level, and 
eventually, depending on the integrated luminosity, these 
measurements could extend out to 8 or even 9 GeV^. 

In order to demonstrate the feasibility of the method 
for realistic cuts, the Monte Carlo event generator PHO- 
KHARA has been extended to simulate this reaction for 
and 71717 final states. Examples for event rates and 
for angular distributions have been presented, which in- 
clude realistic cuts and which demonstrate the feasibility 
of these measurements. NLO corrections to ISR amount 
to typically 1-3%. They are part of the present event gen- 
erator and should be included in a realistic simulation. 

Angular distributions allow the separation of electric 
and magnetic form factors. The general form of these dis- 
tributions has been studied. They become relatively sim- 
ple in the hadronic rest frame, with the z-axis aligned with 
the photon direction. For a specific ("optimal") orienta- 
tion of the 2:-axis the leptonic tensor can even be diagonal- 
ized, which leads to a particularly simple form for the an- 
gular distribution. In fact, this form is quite similar to the 
one observed in the direct electron-positron annihilation 
reaction. The simulation also demonstrates that the sepa- 
ration of the squares of electric and magnetic form factors, 
respectively, can be achieved over a fairly wide Q^-range, 
even if realistic acceptance cuts are imposed. To get access 
to the relative phase between electric and magnetic form 
factors, the determination of the nucleon spin is required, 
which is an attractive option for AA production. We will 
come back to this possibility in a subsequent study. 



with s-y = sin 9-y and c-y — cos 6-y accordingly. We perform 
the following transformations: a rotation of the coordinate 
frame around the a;-axis, such that the new z-axis points 
into the negative photon direction and subsequent boost 
from the laboratory frame along the new z-axis with 7 = 
(s + Q^)/2^y sQ^ into the hadron rest frame. Then, one 
finds 



M V 
Pi, 2 = ^7 




/ 



N 




(22) 



In this coordinate system the space components of the 
leptonic tensor are given by 



ith 



■Via = 




(23) 



(24) 



In combination with the hadronic tensor this leads im- 
mediately to Eq. H1U|) . The leptonic tensor, as given above, 
is evidently symmetric. The combination 



A The diagonal leptonic tensor 

Let us start with the leptonic tensor (see e.g. Eq. (4)) 



in the limit tr^ ^ s and m'^/ s ^ 9: 



Q^ym 

4Q2 



(20) 



where y^^ = ""^ (1 =F cosQ-^). 

Terms proportional to and do not contribute 
as a consequence of current conservation, Qf^Jl^g^^ — 0, 
and have been dropped. The momenta in the laboratory 
frame, with the z-axis pointing along the positron beam 
are given by 





/ 1 







2 










2^ 



I 7^ \ 



V -C7 J 



(21) 



{v\vi+vW2)/'^ = 





4/7' -S7C7/7 I (25) 
-S7C7/7 (/5^ + 4) 



can thus be diagonalized by a rotation of the coordinate 
system around the x-axis, now, however, in the hadron 
rest frame. The rotation angle Ob is obtained from 



tan 9d 
or, alternatively 



A 



A + a 



7(A-a) 



tan(2(?z,) = 

7a 



where 
A = 
with 



X{x,y,z) = + + — 2{xy + xz + yz). 
The eigenvalues of {v\v-[ + v\v2)l2 are given by 
A± = (/3' +d + s^h^ ± A)/2 



(26) 



(27) 



= {P^ + ^-slh^), (28) 



(29) 



(30) 
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and the leptonic tensor in the new "optimal" frame sim- 
phfies to 
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-2— 72diag(0, A_,A^ 



(31) 



Choosing this new "optimal" frame for the definition 
of the baryon direction with angles denoted by 9 and ip, 
only quadratic terms in cos 9 and sin 9 sin (p remain in the 
angular distribution. The combination LijH^^ is now given 

by 



(47ra 



i2 r 



Q^yiy2 



X [ A{\+ + A_) + P%B{X_ sin^ ^sin^ <^ + A+ cos^ 9) 



(32) 

where A and B characterize the hadronic tensor 

H^. = A{Q^Q, - g^,Q^) + ABq^q, . (33) 

For two-particle production this is the most general form 
of the symmetric hadronic tensor, if one sums over polar- 
izations of the particles in the final state. For the nucleon- 
antinucleon final state specifically, one finds 



A = 2\G 



N |2 
Ml 



and B = — 



2t 

r- 1 



G 



N |2 



--IG 



N\2 



.(34) 



Let us emphasize that this choice of coordinates also leads 
to an extremely simple angular distribution for pion pairs, 
and quite generally for any single-particle inclusive distri- 
bution. 

In the the limit of ^ s the rotation angle 9d van- 
ishes, the direction of the photon can be used as (negative) 
z-axis and the leptonic tensor reduces to 



(47ra)2 (1 + 4) 



ym 



and 



(47ra)^ 

Q'^yiy2 



l + cf 



diag(l,l,0) , 



2 A + B(}% sin^ 9 



(35) 



(36) 



which for A and B from Eq. 134|) gives 



Kl^fl + cos^l 



-\G 



N\2 



sm 



(37) 



which is closely related to Eq. 1)12(1 . 
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